Vortices and ring dark solitons in nonlinear amplifying waveguides 
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We consider the generation and propagation of (2+l)-dimensional beams in a nonlinear waveguide 
with the linear gain. Simple self-similar evolution of the beams is achieved at the asymptotic stage, 
if the input beams represent the fundamental mode. On the contrary, if they carry vorticity or 
amplitude nodes (or phase slips), vortex tori and ring dark solitons (RDSs) are generated, featuring 
another type of the self-similar evolution, with an exponentially shrinking vortex core or notch of 
the RDS. Numerical and analytical considerations reveal that these self-similar structures are robust 
entities in amplifying waveguides, being stable against azimuthal perturbations. 
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The topic of self-similarity has been in the focus of in- 
tense research interests in versatile areas of physics p| . In 
nonlinear optics, the intrinsic self-similarity of the non- 
linear Schrodinger (NLS) equation has led to the devel- 
opment of the symmetry-reduction method, which has 
been widely applied to search for exact and asymptotic 
self-similar solutions Such optical waves {similari- 

tons) possess many attractive features that make them 
potentially useful for various applications to fiber-optic 
telecommunications and photonics, as they can main- 
tain the overall shape, while allowing their amplitudes 
and widths to vary, following the modulation of system's 
parameters - the dispersion, nonlinearity, gain, inhomo- 
geneity, and others. Recent studies reveal that similari- 
tons also exist in other fields of ph ysics, i ncluding Bose- 
Einstein condensates and plasmas [10l4l2| . 

Among various types of self-similar modes, there has 
been an increasing interest in the study of asymptoti- 
cally exact parabolic similaritons, since their first experi- 
mental realization in normally dispersive fiber amplifiers 
. One important advantage of these waves is that they 
are free of collapse and filamentation at high power lev- 
els. Another remarkable property of the similaritons is 
that the corresponding output profile is completely de- 
termined by the input power, i.e., all initial profiles with 
the same power evolve toward the same parabolic simi- 
lariton 0,0, which is very attractive for applications 
such as the pulse (or beam) shaping or optical regener- 
ation This conclusion, however, is correct only 
when the initial beams are of the fundamental type, i.e., 
they do not have phase slips or phase singularities, and 



their amplitude profiles do not have nodes. If the in- 
put beam carries a phase slip, a dark soliton in the self- 
similar parabolic background is generated [l6j . On the 
other hand, a spatiotemporal vortex torus will emerge if 
the input beam has an embedded vorticity [T7l - [l9| . 

In this work we study the generation and propagation 
of optical beams inside the two-dimensional nonlinear 
waveguide amplifier, focusing on effects of initial phase 
singularities, phase slips, and amplitude nodes. The non- 
linearity of the medium is assumed to be self-defocusing, 
with refractive index n = no — n2l, where / is the beam's 
intensity, and rt2 is positive (the cubic self-defocusing is 
possible in semiconductor waveguides). The respective 
governing equation for the paraxial optical beam in such 
an amplifying waveguide is the (2-fl)-dimensional NLS 
equation, 

.du 1/9^ 9^ \ , ,o .g 
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where the beam's envelope u, propagation distance z, 
spatial coordinates x and y, and gain g are normalized 
by (fco"2iz5)^^^^, Ljj, Wo, and L~jj^, respectively, with 
ko = 27rrto/A, Ld = koWQ, and wo being, respectively, 
the wavenumber at the input wavelength A, diffraction 
length, and a characteristic transverse width. 

Our first aim is to investigate the self-similar evolution 
of the optical beam in the asymptotic limit of z — > oo by 
considering the intrinsic self-similarity of Eq. ([Ij. For 
the evolution to be exactly self-similar, the functional 
form of the beam's intensity profile must remain un- 
changed at different propagation distances. Accordingly, 
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one can express the beam's intensity as \u{z^x^y)\'^ = 
eMJ^9iz')dz']\U{X,Y)\^/l\ where {X,Y) = {x,y)/e 
with £{z) being a positive function of the propagation 
distance that characterizes the evolution of the beam's 
width, and \ U {X, F) p is a functional form satisfying the 
energy conservation condition, obtained after the sub- 
traction of the gain effect: 



J J \U{X,Y)\'^dXdY 



Uo, 



(2) 



with Uq being the input power of the optical beam. With 
this intensity shape, the expansion velocity of the beam's 
intensity profile is v = {dt/dz)r, where r = {x,y}. 
This vel ocity is equal to the gradient of the beam's phase, 
V — V(/) [20]. Therefore, the phase of the beam undergo- 
ing the self-similar evolution is = (2£)~^ {di/dz)r'^ — 
(/)o(z), where ipoiz) is a phase offset. Thus, by introducing 
the following self-similar transformation. 



u{x,y,z) 



U (X, Y) 



■ exp 
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we cast Eq. ([T]) into the form of 
1 
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where = d^/dX^ + d^/dY^ and =X^ + Y^. 

Without the loss of generality, we let (d-^i/dz^)i^ — 
expljg g{z')dz'], and {d(/)o / dz) i'^ — fiexp[J^ g{z')dz'], 
where /i is a positive constant to be determined. Note 
that the coefficient in front of the diffraction term in Eq. 

exp[— g{z')dz']/2, vanishes at z — ^ oo. There- 
fore, if U {X,Y) is smooth with finite value of V^C/, 
the diffraction term may be neglected in the asymptotic 
regime. The above consideration shows that the exact 
self-similar evolution of the optical beam can be achieved 
in the asymptotic limit, under which Eq. (|4]) yields 



\U\' 
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(5) 



for < 2/i, an d U = otherwise, where the positive 
constant /x = ^JVofn is determined by Eq. Note 
that this profile is essentially the same as produced by 
the Thomas-Fermi approximation for the ground-state 
solution of the two-dimensional Bose-Einstein conden- 
sates described by the Gross-Pitaevskii equation with the 
isotropic parabolic potential and repulsive interatomic in- 
teraction (21)]. Thus, the exact self-similar beams are ro- 
bust entities. 

Hereafter we assume g to be a constant, hence the 
variable characterizing the change of the beam width is 



l{z) 



exp 



(9' 



(6) 



and the phase offset is 



flexp(f) 



(7) 



From Eqs. ((3)) and ([6]), it follows that the amplitude and 
the half width of the exact self-similar beam increase ex- 
ponentially as ,JJIg ex^{gz / A) /2 and \/8ii/g exp{gz/A), 
respectively. As the optical beam expands, a phase-front 
curvature (chirp) develops. Using the stationary-phase 
method [22|, one can find that the spatial spectrum of 
the optical beam is parabolic too. Note that the asymp- 
totically exact self-similar evolution of the optical beam 
is predicated upon the assumption that U is smooth, so 
that the relative strength of the diffraction becomes negli- 
gible when compared to that of the nonlinearity 6] . This 
condition does not hold if the input beam carries vortic- 
ity, or nodes in its amplitude profile, or phase slips. As 
said above, the generation of vortex tori and ring dark 
solitons (RDSs) is expected in that case. 

Now we proceed to the investigation of the dynamics of 
a vortex torus of topological charge S in the asymptotic 
limit. As the local beam's intensity vanishes at the center 
of the vortex, the relative strength of the diffraction in the 
vortical core is much larger than that of the nonlinearity. 
The size of the core range is determined by the coherence 
length, L = S^/l/2\u\^. In the presence of the gain, 
the beam's amplitude increases exponentially, therefore L 
decreases at the same rate. Therefore, at the asymptotic 
stage, the vortex becomes very narrow when compared 
to the whole beam, which makes it reasonable to seek for 
a solution to Eq. Q in the asymptotic limit by setting 



U = ^'(i?)exp(i56'). 



(8) 



where 9 is the azimuthal angle, and neglecting terms con- 
taining gradients of ^I'. Then, we obtain an asymptotic 
expression for the local intensity of the vortex torus. 
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where p = (S'^/2) exp{—gz), and /i in Eq. ([5]) is replaced 
by /ivort such that the phase offset keeps the form of Eq. 
(O with /i replaced by /Uvort- 

From Eq. it follows that, in this approximation, 
the power in the core of the vortex vanishes at R^ < 
/ivort — \/ Mwrt ~ 2p « p//ivort (in the exact solution, the 
beam's intensity vanishes as i?'"^' at i? ^ 0). On the 
other hand, the last term in Eq. © may be neglected 
at large R, hence the beam's intensity profile keeps the 
parabolic form outside of the core. Note that the local 
beam's intensity of the vortex torus at large distances 
is slightly greater than that of the fundamental soliton 
[Eq. dS])] for the same input power, since the optical field 
is removed in the narrow core of the vortex. Therefore, 
the difference of /ivort, which is generated by the local 
nonlinearity, from /i is almost negligible. Thus it can be 
concluded, from the above expressions, that the radius of 
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FIG. 1: (a) The self-similar evolution of the vortex torus of 
topological charge 5 = 1 in the waveguide amplifier, starting 
with input beam u(0,r) = r exp(— exp(i6'), for g = 2. 
From the top to bottom, the propagation distance is 5.25, 
5, 4.75 and 4.5, respectively, (b) the beam's amplitude in 
the core of vortex, which is a linear function of r. (c) The 
phase (phase offset (jjoiz) is ignored) of the vortex torus at 
propagation distance z = 5.25, with the inset showing the 
phase slip tt in the x direction, to confirm that the topological 
charge is 1. Here and in following figures, the solid lines and 
circles represent results of numerical simulations of Eq. ([l]) 
and the analytical prediction, respectively. 



the vortex decreases exponentially, as 



^^5exp(--), 



(10) 



whereas the width and amplitude of the whole beam in- 
crease exponentially, like those of the fundamental spa- 
tial soliton. This analysis predicts peculiarities of the 
self-similar evolution of the vortex torus, in comparison 
with the vortex-free beam. 

It should be emphasized that now U is not only a func- 
tion of X and Y, but also a function of z [Eqs. ^ 
and The latter contradicts to the definition of the 
function U = U{X,Y) in the transformation [Eq. ([3])]. 
This contradiction, however, is negligible: if U is also 
a function of z, Eq. ([4]) will contain an additional term 
ii'^exp[- g{z')dz']Uz. FromEq. jH) it follows that this 

term can be approximated by 2iexp{—gz/2)^^^/^pR~^, 
which is about exp{gz/2) times less than the last term 
in Eq. and hence it is negligible at the asymptotic 
limit. The same argument also holds for the dark solitons 
embedded in the parabolic background shown below. 

Figure 1 displays the evolution of an initial Laguerre- 
Gaussian beam of topological charge = 1 towards the 
vortex torus, whose intensity and phase profiles at the 
asymptotic stage of the evolution are found to be in good 
agreement with the predictions based on Eqs. ([3]), ([HI), dH]) 
and ([5]). Further numerical simulations show that initial 
Laguerre-Gaussian beams with topological charges up to 
S" = 10 evolve into the corresponding vortex tori, indi- 
cating that the latter are robust entities in the waveguide 
amplifiers. 
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FIG. 2: The self-similar evolution of radial dark soli- 
tons in the amplifying waveguide, obtained from in- 
put beams of total power tt with initial profiles (a) 
u{0,r) = M/-^exp(-rV2Vy^)(r2 - W^), and (b) u{0,r) = 
exp(-r^/2H/^) exp{i7r tanh[2(r - W)]/2}. The gain pa- 
rameter is (? = 2 and the initial width is W = 1.5. 



Now, we proceed to the generation and propagation 
of dark solitons in waveguide amplifiers. Since quasi- 
one-dimensional (stripe-shaped) dark solitons eventually 
decay into vortices, we focus on the RDS [23]. Numeri- 
cal simulations reveal that dark solitons of this type can 
be generated from input beams with nodes in the am- 
plitude profile [Fig. 2(a)], or from Gaussian beams with 
ring phase slips [Fig. 2(b)]. In the asymptotic limit, the 
approximate RDS solution of Eq. (j4]) can be written as 



\U\' 



2P 



I - sech^ [^/JlA exp 



(11) 

where A and characterize the depth and location of the 
center of the RDS, respectively. From Eq. pT|) it follows 
that the effective width of the black RDS decreases ex- 
ponentially, as (2/y^) exp(— gz/4), cf. Eq. (fTOl) . The 
character of self-similar evolution of the RDS, which is 
qualitatively similar to that of the vortex, is confirmed by 
simulations displayed in Fig. 2, where one could also find 
that (i) the RDSs generated from the initial phase slip are 
much more shallow than those generated from the initial 
amplitude node, and (ii) the RDSs move outwards due 
to their self-repulsive nature and the expanding parabolic 
background. Finally, Fig. 3 shows the output beam in- 
tensity at z = 5 in the amplifying waveguide with g — 2, 
as obtained from the input beam of power tt, where one 
could see that the vortex and RDSs may coexist, being 
embedded in the self-similar parabolic background. 

Next we investigate the stability of RDSs. By introduc- 
ing the transformation U — ^ J7 exp(— i(/)o), we rewrite Eq. 
dl as iUz + V2[//[2exp(gz)] - \U\'^U = R'^U/2, where 
the new propagation distance Z = {l/2)[exp(gz/2) — 1], 
as obtained from Eq. ([7]). If the coefficient in front 
of the diffraction term is 1 (g = 0), this equation is 
just the governing equation for the dynamics of RDS in 
Bose-Einstein condensates with the harmonic potential 
(23} . and it was demonstrated that (i) in the limit case 
of a quasi-plane soliton, the mass of the deep RDS is 
2 and its oscillation frequency is = -^1/2, and (ii) 
the RDSs persists up to Z ~ 3T = Gtt/lu, and then be- 
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FIG. 3: (a) The coexistence of a vortex with the topological 
charge 5* = 2 and radial dark soliton, (b) the beam's am- 
plitude in the core of vortex (solid lines), where the squares 
represent the numerical fits of |it| as quadratic functions of 



The initial condition is u(0, 



0.309 exp(-r72)(r^ 



3/2)r exp(2i6), and the gain parameter is g = 2. From the 
top to bottom, the propagation distance is 5, 4.75 and 4.5, 
respectively. 



gins to decay due to the azimuthal modulation instabil- 
ity (AMI). When g ^ 0, the effective mass of the deep 
RDS is 2 exp{gz) and its oscillation frequency is approx- 



imately ^/^/2exp{—gz/2). Thus, the RDS would reach 
the AMI threshold at Z' « eV^ir exp{gz/2). In fact, 
because Z is always smaller than Z', we infer that the 
instability threshold will never be reached, in the present 
model, which is indeed conformed by direct numerical 
simulations. Thus, the RDSs are robust entities in the 
waveguide amplifiers, contrary to the common opinion 
that AMI is inevitable for two-dimensional topological 
modes, such as RDS [13,123,123. 

In summary, within the framework of the (2-1-1)- 
dimensional NLS equation including the linear gain, we 
have found three types of self-similar optical beams in 
nonlinear amplifying waveguides, namely, the fundamen- 
tal beams with the parabolic profiles, robust vortex tori 
with different values of the topological charge, and radial 
dark solitons with a parabolic background. These results 
were obtained in the analytical form using the Thomas- 
Fermi type of the asymptotic approximation, and con- 
firmed by numerical simulations. 
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